Scaling of the dynamical properties of the Fermi-Ulam accelerator model 
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The chaotic sea below the lowest energy spanning curve of the complete Fermi-Ulam model (FUM) 
is numerically investigated when the amplitude of oscillation e of the moving wall is small. We use 
scaling analysis near the integrable to non-integrable transition to describe the average energy as 
function of time t and as function of iteration (or collision) number n. If t is employed as independent 
variable, the exponents related to the energy scaling properties of the FUM are different from the 
ones of a well known simplification of this model (SFUM). However, if n is employed as independent 
variable, the exponents are the same for both FUM and SFUM. In the collision number analysis, we 
present analytical arguments supporting that the exponents c* and b* related to the initial velocity 
and to e axe given by c* = —1/2 and b* = —1. We derive also a relation connecting the scaling 
exponents related to the variables time and collision number. Moreover, we show that, differently 
from the SFUM, the average energy in the FUM saturates for long times and we justify the physical 
origins for some differences and similarities observed between the FUM and its simplification. 
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I. INTRODUCTION 

The Fermi accelerator is a dynamical system, pro- 
posed by Enrico Fermi [l| to describe cosmic ray accel- 
eration, in which a charged particle collides with a time 
dependent magnetic field. With base in different appli- 
cations, the original model was later modified and stud- 
ied by other authors. One of them is the well known 
Fermi-Ulam model (FUM) [i,[3|, in which a ball is con- 
fined between a rigid fixed wall and an oscillatory mov- 
ing one. In this model the ball collides elastically with 
the walls and the system is described by an area pre- 
serving map. In high energies regime the phase space 
presents Kolmogorov-Arnold-Moser (KAM) islands sur- 
rounded by locally chaotic regions, which are limited by 
spanning curves. Below the invariant spanning curve of 
lower energy a globally chaotic sea involves KAM islands. 
The presence of invariant spanning curves limits the or- 
bits in phase space impeding unlimited energy growth 
(Fermi acceleration) Another model, very similar 

to FUM, in which the ball is in a gravitational field, so- 
called bouncer 0], presents, differently from FUM, the 
property of Fermi acceleration depending on values of 
control parameter and initial conditions. This difference 
between the two models was later explained by Lichten- 
berg et al Hybrid versions [7|] involving both FUM 
and bouncer and versions of FUM with energy dissipa- 
tion [1] were also explored. Quantum models based on 
FUM and bouncer have also been studied [l^, [Tit . 
The study of such systems is interesting because it allows 
to compare theoretical predictions with experimental re- 
sults Jl2j, JM ' moreover the knowledge about how time- 
dependent perturbations affect the dynamics of Hamil- 
tonian systems is something that needs to be more ex- 



plored. Therefore it is useful to study such perturbations 
in simple models because they furnish insights about 
more complex systems; even more, the formalism used 
in its characterization can be extended to the billiard 
class of problems [H, [H, [M [13 ■ 

The simplified Fermi-Ulam model (SFUM) [2| is an 
approximation in which the position of the moving wall 
is considered as fixed, but it transfers momentum and 
energy to the particle. This geometrical change of the 
complete model can be neglected when the amplitude of 
oscillation is much smaller than the distance between the 
two walls and the velocity of the particle is larger than 
the wall maximal velocity. It is clear that the simpli- 
fied versions can speed up the simulations. However the 
main interest in these simplified models is that they can 
be studied by analytical methods whose results are of- 
ten compared with the numerical results of the simplified 
model. Moreover, these analytical results are, sometimes, 
also useful in study of the full models. 

The moving wall in the Fermi-Ulam model represents 
an external force that acts as a perturbation in the sys- 
tem. If the oscillation amplitude of the moving wall is 
zero the system is integrable but as soon as this ampli- 
tude is different from zero the Fermi-Ulam model behaves 
chaotically [l^. Near the integrable to non-integrable 
transition, average quantities can be described by scaling 
functions. This kind of analysis was originally proposed 
in a study of the SFUM, where the collision number with 
the moving wall is employed as independent variable 19] . 
In a recent work [23], also based on the SFUM, we pro- 
posed a similar analysis in which time is considered as 
independent variable and we showed that the average en- 
ergy can also be described by scaling functions but with 
a different exponents set. Moreover, we showed that the 
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average energy decays at long times. 

In the present work we investigate the chaotic sea be- 
low the spanning curve of lowest energy for the full Fermi- 
Ulam model, where scaling properties of average energies 
are studied on variables time and iteration (or collision) 
number. We show that if the time is employed as inde- 
pendent variable, then the exponents related to the scal- 
ing properties of FUM are not the same that the ones of 
SFUM. However, if the iteration number is employed as 
independent variable, then both FUM and SFUM present 
the same set of exponents. We show that the energy 
decay found for long time in the SFUM [13] does not 
exist for the FUM and we justify physically the origins 
of the similarities and differences between the FUM and 
SFUM. We provide also some analytical results for the 
scaling exponents and, moreover, we show that although 
the exponents related to the variables time and collision 
number are not the same a relation between them can be 
established. 

This paper is organized as follows: In the next sec- 
tion the Fermi-Ulam model and the average quantities 
of interest are defined. We present also the procedures 
to obtain these averages when time is the independent 
variable. In Section III we present the results of the nu- 
merical simulations, the scaling properties of average en- 
ergies as function of both time and collision number, and 
we determine the exponents related to these averages. In 
Section IV analytical arguments to determine some ex- 
ponents and to derive a relation between the exponents 
of time and collision number analyses are discussed. Fi- 
nally, we draw the conclusions in Section V, where we 
also present a summary of this present work. 

II. THE FERMI-ULAM MODEL 

The Fermi-Ulam model describes the motion of a clas- 
sical particle bouncing between two rigid walls, one of 
which is fixed at position x — and other that is 
moving periodically in time whose position is given by 
Xw{t") = xo + e' cos{ujt" + (j)o). Here xq is the equi- 
librium position, e' is the oscillation amplitude, w is a 
frequency, t" is time and (j)o is the initial phase. In order 
to work with dimensionless variables we perform scale 
changes in length = Xw/xq and in time t' — cut". 
With these new variables the system has just one param- 
eter, namely e = s'/xq, and we write the position of the 
moving wall as = 1-1- ecos{t' + 0o). The parti- 

cle moves freely between the walls and collides elastically 
with them. In this manner the FUM can be described 
by a map T{Vn,4>n) = (Ki-i-i, '/'n+i) which gives the ve- 
locity of the particle and the phase of the moving wall 
immediately after each collision [l^ 

T ^ I ^ ^ 2esin((?;>„+i) 

[ (j>n+i = (i>n + Ai„+i mod 2tt ' 

Here the term 2e sin((/)„+i) gives the fraction of velocity 
gained or lost in collision and At„_|_i = — t„ is the 



time between two collisions with the moving wall, which 
is given by the smallest solution of 

KAt„+i-(l+ecos</>„) = ±[l + ecos(Ai„+i-|-0„)]. (2) 

The plus sign in the above equations corresponds to 
the situation in which the particle collides with the fixed 
wall before hitting the moving one {indirect collisions)] 
the minus sign corresponds to the situation in which the 
particle hits successively with the moving wall [direct col- 
lisions or successive collisions). 

Let us define V'^{t') as the square velocity at time t' . 
We are interested in the scaling properties of the dimen- 
sionless energy E = 2 Energy/mcj^XQ averaged over an 
ensemble of M samples that belong to the chaotic sea. 
Such samples are characterized by initial phases 4>q of 
the moving wall randomly chosen in an interval /. If 
the initial velocity Vq is small enough then we can use 
/ — [0, 27r). Moreover we consider that the particle starts 
at the fixed wall position with velocity Vq > 0. Thus, the 
time of the first collision of the particle with the moving 
wall is given by Ti « I/Vq and we define a variable t as 
t = t' — I/Vq. In this way the time t starts at the first 
collision instant and we write the average energy as 

M 

E{t,s,Vo)^-J2v^'{t) , (3) 

where j refers to a sample. 

We consider also another kind of average where the 
square velocity is firstly averaged over the orbit of a sam- 
ple as 

^(0 = ^£ V'ir)dr , (4) 

and then we perform the average of the energy in an 
ensemble of M samples as defined below 

1 *^ _ 

E{t,s,Vo)^-Y.^^Ht) . (5) 

As the particle moves freely between the walls, the 
square of its velocity is constant between two impacts 
with the moving wall and the integral in Eq. Q can be 
numerically evaluated without difficulties. 

We will now describe some numerical procedures used 
to evaluate the averages quantities on time. The dynam- 
ics of FUM as described in Eq. ([T]) evolves in a discrete 
variable, namely the collision number n. As we are inter- 
ested in the evolution of the average energies as defined 
in Eqs. ([3]) and ([5]), where time i is a continuous vari- 
able, we can use the map given by Eq. ([T]) to speed up 
the calculation process. We evaluate the averages in Eqs. 
^ and ^ at discrete, logarithmic spaced, values of time 
t = ti, t2, ■ • ■ , t]sf- We known that at i = 0, the first 
collision instant, the energy of the particle is and, as 
the particle is at the moving wall position, the next col- 
lision takes place at tc = At2, where At2 is obtained by 
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solving Eq. ([2]). To evaluate the energy Ei of the first 
sample of the ensemble at time t = ti we follow the pro- 
cedure: (1) If ti < tc, the energy of the particle at ti will 
be El — Vi- (2) Otherwise ti > tc and a collision occurs 
at time t — tc = At2. Then we employ Eq. ([1]) updating 
the velocity of the particle to V2 and the next collision 
time is given by tc = Ai2 + Ata. 

It means that if case 1 is satisfied then Ei{ti) was de- 
termined. Otherwise we update the next collision instant 
and the velocity of the particle. If now ti < tc (case 1) 
then El = but if ti is still larger than t^ we repeat 
case 2 updating both V and tc- The reasoning is ba- 
sically to repeat the above procedure until ti < tc and 
then to update the energy Ei{ti). We follow a similar 
proceeding, with appropriate time intervals, to evaluate 
£^1(^2), £^1(^3), Ei{tN). 

The procedure is the same for all samples of the en- 
semble giving E2, . . . , Em from ti until t^. Then the 
average energy in Eq. ^ is performed doing E{ti) = 
[Ei{ti) + E2{ti) + ... + EM{ti)]/M, E{tN) = 

[Ei{tN) + E2{tN) + ■■■+ EM{tN)]/M. The average en- 
ergy E{t, e, Vq) in Eq. ^ can be evaluated by a similar 
procedure. 

Since the map in Eq. ([1]) gives the velocity and phase 
after each collision with the moving wall, then the av- 
erages on variable n can be defined in a more direct 
manner. We first consider the average velocity over 
the orbit generated from an initial phase (j)o, defined as 
V'^{n) = ^^-i^Y X]r=o ^i^- Considering an ensemble with M 
samples, characterized by initial conditions that belong 
to the chaotic sea, we define the average energy as 

1 — 

<E>in,e,Vo)^—J2^^'^n) . (6) 

III. RESULTS 

In Fig. [Ija) and [IJb) we show, respectively, the nu- 
merical results for the average energy E{t, e, Vq) when 
the initial velocity is small, or Vq << e, and the situa- 
tion in which Vq » e. The averages defined in Eqs. ([3]) 
and ([5]) were performed in an ensemble with M — 2x10'^ 
samples. Although Vb >> e, the energy curves shown in 
[TJb) are obtained from orbits that belong to the chaotic 
sea. In Fig. [Tfb) we used the same values of e as those 
shown in Fig. [ija). 

As we can see more clearly in Fig. [Ua) the energy is 
constant until a time ti, grows up to a time t2 and then 
reaches a stationary value for large values of time. As 
show in [l^l the energy in the simplified FUM presents a 
slow decay in time. Since in SFUM the oscillating wall is 
considered fixed at position a; = I, the time between two 
collisions with the oscillating wall is 2/y„, and successive 
collisions do not occur. Therefore, if after a collision 
the particle has very low velocity, then it remains for a 
long time, 2/Vm with low energy. In this way, at time 
t many realizations are in such condition originating the 
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FIG. 1: (Color online) Log-log plots of the energy E{t) defined 
in Eq. Q, averaged over M = 2 x 10^ samples, as function 
of time for (a) three values of e and Vo = 1 x 10~® (Vo << e) 
and (b) three values of e, the same ones as those shown in 
(a), and three values of initial velocity Vb >> e. 

slow decay in the average energy ior t » t2. In FUM 
this situation does not occur. If the particle, after a 
collision, losses almost all its energy, a successive collision 
occurs increasing the energy of the particle. Therefore, 
in FUM the average energy does not decay, but presents 
a saturation regime for t » t2. Moreover, the energy 
for t»t2 can be described as 

E{t,e,VQ) oc g{e) , 

9{e) (X . (7) 

The value of exponent (3 can be determined by search- 
ing for the best collapse of the energy curves in the 
asymptotic regime. We performed simulations for val- 
ues of £ between 1 x 10~^ and I x 10^'^ and we obtained 
the average value (3 — 1. 01 ± 0.01. The crossover time t2 
obeys the relation 

i2 oc , (8) 

and the better fit in a plot of t2 as a function of e gives 
z = -1.50 ± 0.03, as shown in Fig. H 

The time <i is an average time in which the second in- 
direct collisions happen. After that new collisions occur 
and the energy i?(i,e, Vb) begins to increase. Moreover, 
ti is related to an initial transient and it affects the be- 
havior of i?(t, e,Vb) for t « t2- Therefore, between ti 
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FIG. 2: Log-log plot of the crossover time as function of 
parameter e. 

and t2 we have a crossover region in which the growth 
exponent can not be directly evaluated. 

For large values of time {t >> ti), the energy 
E{t, e, Vq) can be written as a scaling function, namely 

E{t,e,Vo)^lE{rt,l''e,l''Vo) . (9) 

Here a, b and c are scaling exponents and I is the scaling 
factor. In the limit Vb << e (Fig- [Ha-)) we can choose 
I = e~^l^ and write the above equation as E{t,e,0) = 
e-i/''£;(e-°/4, 1,0) cx £-'^/'' fie-^/H). For t » t2 this 
relation can be written as E{t, e, 0) cx e^^^''. Then, from 
Eq. ([7]) we have that (3 = —1/6. Using the simula- 
tion value P = 1.01 ± 0.01 we obtain b = -0.99 ± 0.01. 
From Eq. ([8|) we derive the relation z = a/b. Since 
the simulations furnish z ~ —1.50 ± 0.03, it follows that 
a ^ 1.49 ± 0.04. 

Fig. [3]Ja) shows the rescaled energy E{t,e,Vo)/l as 
function of rescaled time We can see that with these 
new coordinates the energy curves, originally depicted in 
Fig. [ija), collapse onto a universal curve in the limit of 
long time, after an initial transient. We emphasize that 
the scaling behavior is valid only for small values of e. 

The argument to obtain the c exponent, namely c = 
&/2, is presented in the next section. Using this relation 
we obtain that c = -0.495 ± 0.005. Fig. ^h) shows 
the collapse of the curves depicted in Fig. [Hb) for the 
appropriate chosen initial velocities. 

Depending on initial velocity and phase, the velocity 
after a collision with the moving wall can be very small in 
such way that a direct collision occurs and gives an imme- 
diate increase in energy. This scenario is more common 
for small values of Vq and t and does not allow us to 
find a scaling description to the initial transient at small 
values of t. However, this transient is important because 
it affects the power-law growth of the energy between ti 
and ^2- 

We can use Eqs. ([7]) and ^ to describe the average en- 
ergy E{t,e,Vo), defined in Eq. (O, as a scahng function 
just changing E by E. Figs. Ua) and[Ub) present the en- 
ergy E{t,e,Vo) as function of time i for Vq « £ (small 
initial velocity) and Vq >> e, respectively. The values 
of e are the same as those in Fig. [Tfa). In Figs. O^a) 
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FIG. 3: (Color online) The rescaled energy E{t)/l as function 
of rescaled time tl"' for three values of e in limit of long time. 
In (a) it is shown the collapse for Vo << £ (Vo = 1 x 10~®) and 
in (b) the collapse for three values of initial velocity Vo >> e. 
We chose / = e~^^^ and the exponents are a = 1.49 ± 0.04, 
h = -0.99 ± 0.01 and c = -0.495 ± 0.005. 

and [5[b) it is shown the rescaled energy as function of 
rescaled time in limit of large t {t >> ti). Following the 
same reasoning employed in analysis of the average en- 
ergy E{t, e, Vo), we obtain the exponents a — 1.50±0.04, 
b = -0.995 ± 0.003 and c = -0.498 ± 0.001. Consid- 
ering the uncertainties we observe that both averages 
E{t,e,Vo) and E{t,e,Vo) are described by scahng func- 
tions with basically the same exponents set. Therefore, 
we will now use the average exponents a — 1.50 ± 0.04, 
b = -0.993 ± 0.007 and c = -0.497 ± 0.003. Note that 
these values are different than the ones of the SFUM 
[13, namely, az = 1.35 ± 0.05, 63 = -0.90 ± 0.03 and 
C2 = -0.45 ±0.01. 

Fig. [ni^a) shows the average energy < E > {n,e,Vo) for 
two values of e and different initial velocities, including 
the situations (i) Vb << e and (ii) Vq » e. Now we 
have that 

<£;> (n,e,I/o) = > (r%i, Vo) ■ (10) 

We can determine the exponents a* and b* from the en- 
ergy curves for Vq << e. When Vq — 1x10^^, we observe 
in Fig. [S]Ja) that the average energy < E > {n, s, Vq) 
presents two regimes. For n < the energy has a power- 
law growth and for n >> Ux the energy is constant. 
Moreover, we obtain that oc with z = —1.01 ±0.02. 
Since in the limit of large n the energy depends only 
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FIG. 4: (Color online) The average energy E{t), defined in 
Eq. (O, as function of time t for three values of e and (a) 
three values of initial velocity Vo << £ and (b) three values 
of velocity Vb >> e. The averages were performed with M — 
2 X 10^ samples. 



FIG. 5: (Color online) Log-log plots of the rescaled energy 
E{t)/l as function of rescaled time tl"' in limit of large time for 
three values of e. We have that (a) Vb << £ and (b) Vb >> e. 
We chose I — e~^^^ and the exponents are a = 1.50 ± 0.04, 
h = -0.995 ± 0.003 and c = -0.498 ± 0.001. 



on e, we have that < E > (n,e,Vb) = g{e) ex with 
(3 = 1.02 ± 0.02. We can choose / = e''^/''' and rewrite 
the above equation as < i? > (rt, e, Vb) = e^^^'' < 
E > (e~° n, l,e^'^'/'''Vb). From this relation we ob- 
tain that z = a* /b* = —1.01 ± 0.02. In hmit n >> 
we can write < E > {n, e, Vb) oc e^^/^ . Therefore, we 
have that (3 = -1/5* = 1.02 ± 0.02. This implies that 
b* = -0.98±0.02 and a* = 0.99±0.04. We use a connec- 
tion between the simplified FUM and the standard map, 
described in the next section, to obtain the value of the 
exponent c*, namely, c* = 5*/2 = -0.49 ± 0.01. 

When Vb >> £ we observe in Fig. [DJa) that the energy 
curves present two characteristic iteration numbers, 
and n". We can also observe that rt" « for Vb << e. 
Therefore, we must consider two situations: n" << n^, 
for small initial velocities (Vb << e), and n" ~ n^, for 
Vb >> e. The energy curves with Vb = 1 x 10^'' (n" w 0) 
present only two regimes: (1) a power-law growth for 
n « n'^ and (2) a saturation regime for n » n'^. On 
the other hand, for the energy curves with Vb = 1 x 10""^ 
and Vb = 10^/2 X 10-3 « < n'J shown in Fig. Eta), we 
have three regimes: (1) the energy is basically constant 
for n << n", (2) for n" < n < the energy grows and 
begins to follow the curve of Vb = 1 x 10~^ and (3) the 
energy curves reach a saturation regime for n >> n'^. 

In Fig. IHb) it is shown the rescaled energy < E > /I 
as function of the rescaled interactions number nl"" . As 



TABLE I: Numerical estimations of the scaling exponents for 
both FUM and SFUM. The exponents a, b and c describe the 
scaling properties of the average energy as function of time t 
while the scaling relations of the average energy as function 
of collision number n are characterized by the exponents a* , 
b* and c*. The uncertainties are given between parenthesis. 





FUM 


SFUM 




t 


n 


t \20} 


n \19] 


a, a* 


1.50(4) 


0.99(4) 


1.35(5) 


0.99(3) 


6, 6* 


-0.993(7) 


-0.98(2) 


-0.90(3) 


-0.977(6) 


c, c* 


-0.497(3) 


-0.49(1) 


-0.45(1) 


-0.489(3) 



we can observe, the energy curves, after a small initial 
transient, collapse onto a universal curve, even for Vb >> 
e, with the exponents a* = 0.99 ±0.04, b* = -0.98 ±0.02 
and c* ~ —0.49 ± 0.01. It is important to note that this 
set of exponents is the same, within the uncertainties, to 
that of the SFUM [l^. The exponents of the average 
energies as function of both time t and collision number 
n for the FUM and its simplification are shown in Table 

in 
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FIG. 6: (Color online) (a) Average energy for two values of e 
and three values of initial velocity, Vo, as function of collision 
number, n. (b) The rescaled average energy as function of 
rescaled collision number shows that, after an initial transient, 
the energy curves collapse onto a universal curve. 



and (iii) a set of terms ±£^ < sin((/)i) sin(0j) >. Observe 
that < sin(0j) > and < sin(0i) sin(^j) > take values 
in the interval [—1, 1]. Since we are in the region below 
the first spanning curve, the maximal initial value Vb,max 

Therefore, for small e, 



must be of order of V* 



U/2 



^0 



< 



max ^ ^ 

< > 



e and n small enough we have that 



> 



^0 ' 



implying that < >f 



^0 ■ 



Assuming that the scaling relation is valid in this limit, 
we obtain the relation < >« 1^+^'^' < Vj^ >, which 
furnishes c* = ~l/2. It is worth mentioning that the 
numerical results for Vq >> £ shows that the scaling 
occurs in the beginning of the simulation without any 
transient (see Fig. ([6]). Moreover, using the result c* = 
b*/2, we also obtain that b* = —1. 

Finally we present an argument that relates the ex- 
ponents a* and a, the scaling dimensions of variables n 
and t, respectively. As the average energy presents a 
saturation regime at large values of time, we can define 
the average velocity V = ^J~E and Ai as the average 
time between two arbitrary collisions. We can also de- 
fine L as the average distance that the particle travels 
between these collisions in such way that we can write 
V iSAjL = An, where An is the average collisions number 
with the moving wall that take place in the time inter- 
val Ai. In terms of the rescaled variables we have that 
F'Ai'/L' = An', or l^Vl^M/L = Z"^' An. AsL'^L we 
have, therefore, that the exponent of collision number, 
a*, is related to the exponent of time, a, hy a* = a + c. 
This result is in good agreement with the simulations. 



IV. ANALYTICAL ARGUMENTS 

Let us first derive a relation between the exponents c 
(c*) and b (6*). We follow the same lines of Leonel et al. 
[T8| . Performing the variable change /„ — 2/V* + 2{V* — 
Vn/V*"^), where V* is a typical velocity near to the lowest 
spanning curve, and a linearization around V* , the FUM 
transforms into a standard map which is described by 
In+i = In- Kcs sin 0„+i and (/)„4.i = 0„ -I- /„. Here K^s 
is an effective control parameter given by K^s — 4e/F*^. 
Note that the Standard model presents a transition from 
local to globally stochastic behavior a,t K — Kc ^ 0.972. 
The values of V* in the lowest spanning curve of the 
FUM furnish a Kcs with value approximately the same as 
Kc, independent of e. Then, we use the scaled variables 
e' = I'^e and V*' = l'=V* to obtain K^b = 4e'/y*'^ = 
4(Z''e)/(ry*)2. This implies that c = 5/2. 

Now we present an heuristic argument to support that 
c* = —1/2. From Eq. ^ we can write, for n = 1, that 

= ± 4eVbsin(0i) + Ae^ sin^ {(Iji) . (11) 

A similar equation can be obtained for V2 by changing 
Vo by Vi and by 02- Then, we replace Vi by Eq. 
pT|) . This iteration procedure can be done for with 
arbitrary n. If now we take the average in the ensemble of 
initial phases < V^ >, we always find a sum of three kinds 
of terms: (i) Vq, (ii) a set of terms ±Voe < sin(0j) > 



V. CONCLUSIONS 

We studied the scaling properties of the chaotic sea be- 
low the lowest energy spanning curve of FUM considering 
average energies as function of time t and collision num- 
ber n. In limit of large t {t >> the average energies 
E and E of FUM can be described by scaling functions 
with exponents a « 3/2, 6 « — 1 and c « —1/2 (see Table 
HI). This values are different than the ones of the SFUM 
namely, aa = 1.35 ± 0.05, 62 = -0.90 ± 0.03 and 
C2 = —0.45 ± 0.01. The scaling descriptions of the aver- 
age energies are also hold on variable n, with exponents 
a* « 1, 6* « —1 and c* « —1/2. These values are basi- 
cally the same as those of the SFUM [l9| . We employed 
some analytical arguments to determine the exponents 
c* — —1/2 and b* = —1. We observe also that, in the 
full model, the scaling exponents related to variables e 
and Vo are, within the uncertainties, the same for both 
n and t analyses, or 6* « 6 and c* « c. We can also note 
that the exponents related to the variables t and n (a and 
a* , respectively) are note the same. However, we provide 
an heuristic argument that establishes a connection be- 
tween these exponents by the relation a* — a + c. It is 
important to stress that the scaling analysis only holds 
for small values of the rescaled amplitude e, close enough 
to the integrable (e = 0) to non-integrable (e 7^ 0) tran- 
sition. Moreover, we observe that the energy of the FUM 
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is constant up to a time ti, grows to a time t2 and, then, 
reaches a stationary value. In the SFUM the average en- 
ergy, differently^ of FUM, presents a slow decay for large 
t {t >> p3]. This striking result is a consequence 
of the approximation used to define the SFUM. As the 
oscillating wall is considered fixed in space, the time be- 
tween collisions is 2/V and successive collisions do not 
occur. Therefore, eventually after a collision with the 
moving wall, the particle has very low velocity V and re- 
mains for a long time (2/y) with low energy, originating 
a slow decay in energy for i >> <2- In the full model 
it is different: if after a collision with the moving wall 
the particle losses almost all its energy, then a successive 
collision occurs increasing the energy of the particle and 
the decay in energy at t >> t2 is not observed. The 
analyses on variable n furnish basically the same results 
for both FUM and SFUM because between one arbitrary 
collision and the next one we always have that Art = 1, 
independently if the collisions are direct (successive) or 
indirect. 

Summarizing, we employed scaling analyses to describe 
the properties of average energies as function of t and as 
function of n in regime of small amplitudes of oscillation 
of the moving wall, e « 0, for the full Fermi-Ulam model. 
We observe that i) the scaling exponents related to the 
variables e and Vq are basically the same for both t and 
n analyses, namely, 6 ~ 6* ^ — 1 and c ^ c* « — 1/2. ii) 



We also observe that the scaling exponents related to the 
variables t and n are not the same, 0^3/2 and a* ~ 1, 
respectively, and, by a simple analytical reasoning, we 
show that these exponents are connected by the relation 
a* = a + c. iii) Performing some analytical analyses on 
variable n we derive that h* — —1 and c* ^1/2. iv) 
Considering also the results of the simplified Fermi-Ulam 
model, ^] and [23| , we observe that the scaling analyses 
of FUM and SFUM on variable t are characterized by 
different exponents sets, while the analyses on variable n 
furnish, basically, the same exponents set for both FUM 
and SFUM. v) We observe also that the successive colli- 
sions play an important rule when the Fermi-Ulam model 
is studied on variable t by, mainly, preventing the energy 
decay for long times, vi) In the analysis on variable n di- 
rect (successive) and indirect collisions play a similar rule 
and, therefore, the energy curves as function of n present, 
for both FUM and SFUM, the same behavior and they 
are described by a single set of scaling exponents. 
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